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' Abstract 

^ ■ The effect of quantum particle production near the big rip singularity 

' has been investigated previously, with the conclusion that the energy of 

OO . the produced particle decreases as the future singularity is approached. 

0^ ' Hence, the effect of particle production would not be effective to avoid 

I the big rip singularity. That calculation was performed by introducing 

■ an ultra-violet cut-off. In the present work we consider a renormalization 

' of the energy-momentum tensor, obtainning a different expression for the 

, particle production. The new expression seems to indicate that the effect 

' of particle production may be dominant as the singularity is approached. 

> : 

^ . PACS: 98.80.-k, 04.62. -hv 



1 Introduction 

The particle creation in curved space-time is a quantum effect associated with 
the non-unicity of the vacuum state when the space-time does not have a time- 
like killing vector [l][2]. In this situation the curvature changes with time, and 
the characterization of the vacuum state depends on time. Hence, if we fix an 
initial vacuum state at, let us say, t = 0, this vaccum state becomes later a 
multiparticle state. This is just another way of saying the particles are created 
as the space-time evolves with time. 

There are two main situations where the particle creation as a gravitational 
phenomena has been especially studied: the space-time of a black hole and the 
expanding universe. If we consider the fact that the black hole must be formed 
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dynamically by the gravitational collapse process, there is no time invariance of 
the configuration, the curvature grows with time, and particle must be created 
as the event horizon is formed. On the other hand, the time invariance does not 
exist when the universe expands, and we must also expect that the phenomena 
of particle creation must appear in the usual cosmological scenarios. We must 
expect that the particle creation rate must be linked with the curvature of the 
four-dimensional space-time. 

Considering the usual Friedmann universe, with homogeneity and isotropy, 
there are two situations where the particle creation phenomena does not occur. 
The first one, is when the universe is in a de Sitter state. The de Sitter space- 
time is maximally symmetric. Hence, it contains a timelike Killing vector, and 
the vacuum state must be unique. The other situation is a radiative universe. 
Since the electromagnetic field is conformal invariant, the particle production is 
zero; in fact, an isotropic conformal field does not lead to any particle creation. 

For other types of Friedmann universes, the situation is more complex. In the 
cases where the strong energy condition p -I- 3p > (p the energy density and p 
the pressure) is satisfied, the curvature is initially infinite, decreasing later with 
time. In such a situation, the greastet problem is to define an initial vacuum 
state: such definition becomes quite arbitrary, and the vacuum state appears 
naturally only in the asymptotic future. When the strong energy condition 
is violated, however, it is possible to define uniquely a vacuum state, and the 
problem of particle creation seems to be well posed. 

The problem of particle creation in a Friedmann universe has been analysed 
recently in reference [3]. The main motivation of that analysis is the fact that 
the universe today seems to be dominated by a phantom fiuid [5, a fiuid that 
violates the dominant energy condition p + p > Q. A universe dominated by a 
phantom fluid will end inevitably in a future singurality, in a finite proper time, 
called big-rip: as the universe expands, the energy density of the phantom fluid 
increases, becoming divergent at a given moment in the future [5 . If this is 
so, the question it was tried to be answered in that work was: can the particle 
creation to be so effective that it leads to the avoidance of the big rip? It was 
used a massless scalar field leaving in a universe dominated by a fiuid with an 
equation of state p ^ ap. In particular, it was payed special attention to the 
case a < —1/3, the situation where the strong energy condition is violated. 

The analysis of this problem becomes cumbersome due to the fact that the 
energy-momentum tensor for the scalar field becomes divergent in the ultraviolet 
limit. Hence, a renormalization procedure must be employed. In the case of 
curved space-time this is not a simple task. In reference |3] it was chosen to 
introduce an ultraviolet cut-off connected with the Planckian frequency. The 
reason for that is the fact that such problem involves gravitation, and above the 
Planck frequency we enter in the quantum gravity domain, and very probable 
the dispersion relation connecting the frequency uj and the wavenumber k must 
be changed. 

The answer for the question of the avoidance of the big-rip due to particle 
creation found in [3] using the procedure described above was negative: the 
number of particle created goes to infinity, but the associated energy density 
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goes to zero, as the big-rip is approached. This seems to agree with the result 
obtaining in reference [6], where a similar study was performed in the case 
of the sudden singularity [7j, a singularity that occurs also in the future but 
without violation neither of the strong energy condition nor the dominant energy 
condition: the sudden singularity is robust with respect to the particle creation 
phenomena. 

However, it must be remarked that, in contrast with what happens in the case 
of reference [5] , in the sudden singularity scenario the subtraction of the infinities 
in the energy-momentum tensor is a much easier task, and is possible to obtain a 
renormalized energy- momentum tensor (ignoring all transplanckian problems). 
Here, we would like to come back to the problem treated in reference [3], and 
ignore the transplanckian problem also in big-rip scenario, trying to renormalize 
the energy-momentum tensor for the massless scalar field. We employ the n- 
wave method, which is quite convinient for an isotropic universe. However, this 
task becomes more complicated due to the existence of logarithmic divergence 
in the ultraviolet limit. We propose a simple way to deal with this problem, 
without introducing any new arbitray parameter. We obtain a finite energy- 
momentum tensor. However, using now this renormalized energy, we verify that 
the energy associated with the particle production diverges near the big rip, and 
it can become the dominant energy component, altering the final state of the 
universe. The unicity of the answer and its modification due to transplanckian 
considerations remain open questions. However, our result points out that the 
problem can be much more complex than it seems at first glance. 

This paper is organized as follows. In next section we revise the general 
problem of renormalization of a quantum scalar field living in a FRW space- 
time. In section 3 we apply the formalism to the case of a radiative universe, 
and in section 4 to the case of a de Sitter universe. In section 5 the general 
inflationary case is discussed, showing that the method has some limitations 
due to the appearence of an ultraviolet divergence. A full renormalisation is 
proposed in section 6, obtainning as result that the big rip may be avoid due to 
quantum effects. In section 7 we present our conclusions. 

2 General problem of renormalization 

The Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 



describes a spatially flat, expanding universe, with a{t) being the time-dependent 
scale factor. In terms of the conformal time defined by dt ~ adi], the FLRW 
metric takes the form 



(1) 



ds^ = a^{dT]^ ~ dx^ - dy^ - dz^). 



(2) 
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If the universe is filled with a barotropic perfect fluid, with the equation of state 
p = ap, the Friedmand equation 



a'V 8ttG 



pa' (3) 



a J 3 

(primes meaning derivative with respect to the conformal time ?/) implies 

a^V^, P^Y^- (4) 

For a > —1/3, the strong energy condition p + 3p > is satisfied and the 
universe expands desacelerating as r/ — > oo; for a < —1/3 the universe expands 
accelerating as ry — > 0_. In particular, a = —1 corresponds to a de Sitter 
universe, while a = 1/3 implies a radiation-dominated universe. If a < —1, the 
fluid is called phantom, and the universe must faces a big-rip in the future. 

The fundamental equation for a massless scalar field is living in a FLRW 
space-time is, 

cl>'l + + k^<j>k = 0. (5) 

a 

The general expressions for the energy and pressure of the quantum scalar field 
are given by the following expressions: 

P^^ l^°"{^^'k<P'l + k'^krkY^dk; (6) 

These expressions can be written as 

P = I k'Ekdk, i?fe = l|0',<^'*+fc2^fe0*,|, (8) 
P = k^Pkdk, Pfc = -|0',.<^'*-_0fe0,..j. (9) 

In general, the above expressions for the energy and the pressure have quar- 
tic, quadratic and logarithmic divergencies in the ultraviolet limit, while it is 
regular in the infrared limit. In order to cope with these divergencies, we em- 
ploy the n-wave method described in the reference [S] . Essentially, this method 
consists in subtracting terms obtainned by expanding E). and Pk in powers of 

Er^Ek-E^k-El-^El (10) 

pr^Pk-p^-p^-lpl (11) 

where 

El = hm J , Pi - hm " \ , (12) 
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with the definitions, 



E]: = -Ekink), Pfc" = -Pkink). 
n n 

It will come out more convient to express the derivatives as 

df 



dr 



dn df _ _rf_dl 

dnr^ dn 2 9n ' 



(13) 



(14) 



and subsequently for the higher derivatives, / being either E^ or P/j. Hence, we 
find: 



df 

d^f 



~Ydn' 

4 dn 4 dn^ 



(15) 
(16) 



It will be necessary later to use the following expressions for the Hankel's 
functions in the limit of large values for the argument: 



Hi^\x) ~ ^i^P[u,x)-iQ{v,x)^e 



(17) 
(18) 



where x 



v + 



and 



fe=o 



(2a; 



\2k 



2\{8x) 



+ 



(4i/2 - l)(4!/2 _ 9)(4zy2 _ 25)(4i/2 - 49) 



Qiv, x) 



E(-i) 



4!(8a;)4 

{v, 2k + 1) 4!/2 - 1 



+ 



(19) 



(2.t)2'=+i 



8a; 



3!(8a;)3 
r(l/2 + n + fc) 



+ 



fc!r(l/2 + n-fc)' 



(20) 
(21) 



Wc will consider from now on two particular cases, the radiative case and 
the de Sitter well as the general inflationary and phantom cases. 
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3 Radiative case 



Let us consider now the particular case of the the flat universe dominated by the 
radiative fluid. Since the radiative fluid is conformal invariant, and the universe 
is isotropic and homogenous, we must expect that the rate of particle creation 
is zero, as it is stated in reference [S] (see also reference [S])- 

For a radiative universe, the scale factor behaves as a = aorj. Hence, the 
equation ([5]) reads now. 



2— + fc^0 = 0. (22) 
V 



The general solution is 



A[kr} — k-x) A(kri-\-k-x) 

c^^ci{k) +C2(fc) . (23) 



We find the initial Bunch-Davies vacuum if, for example, 

cr = /^, (24) 
Now the energy density and the pressure read, 

» - ^f{2>^'.^}^^k-^j;{2>.'.^}^.^^, (25) 



These expressions for the energy and pressure obey the conservation law. 



p' + 3-(p + p)=0. (27) 



2. . 11 1 



From (EEl), we find: 



There is a quartic and a quadratic divergences both for the energy and for the 
pressure. 

In order to give sense to the energy and pressure expression, let us proceed 
with the renormalization schema described in section 2. Using the expressions 
for e P^, we obtain the following relations: 

^" = 2^ El = ^, El^O, (30) 
P'.-~.^ Pl-j^^ Pl = ^ (31) 
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Hence, we have, 

Sr = Ek-El-El-^El^Q, (32) 
Pr = P,-P'^-Pl-\PI=Q. (33) 

Consequently, the initial expression for the energy diverges but, after renor- 
malization, it is equal to zero, in agreement with the fact that no particle is 
created during the radiative phase. 

Note that the equation (P^ . after the redefinition = f , takes the form 

r + = 0. (34) 

From this expression it is clear that we could expect no particle creation from 
the begining. 



4 Inflation: de Sitter phase 

In a de Sitter phase, resulting from the equation of state p = —p {a = —1), the 
scale factor is given by 

a = --. (35) 
1] 



The Klein- Gordon equation reads. 



2^ + k^(f) = 0. (36) 



The solution is, 



^ = ??'/'|ci H^%{kv) + C2 H^%{kv)j. (37) 

Notice that, from now on, wc have made the substitution rj —>■ — in order to 
avoid the rcpctivie use of the absolute value of the original rj parameter which 
is defined in the negative real axis. 

The initial condition is imposed at — s- cx), where 



We can obtain an initial Bunch-Davies vacuum state if 

ci = C2 = 0. (39) 

Hence, we have the final solution 

= 773/2ciijW(fc^). (40) 
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Using the recurrence relations for the Hankel's functions, we obtain the 
following expression for the energy and the pressure: 

p = 4ncW l\''I^H[%{krj)H[%{kr]) + H^)l{kv)Hl^^^^^ (41) 
p = AncW fc'{i?l;^(M Hfj>^{kii) - ii?i;^(M iff/) (fcr?)|dfc. (42) 

These expressions contain a quartic and a quadratic divergences as in the radia- 
tive case. In fact, these expressions can be simplified using the following forms 
for the Hankel's functions: 

<JW = -<\/J'" . f = (43) 



With these expressions for the Hankel's functions, the expression for the energy 
and the pressure reduces to: 

We see to appear the usual quartic and quadratic divergencies in the ultraviolet 
limit. There is no divergence in infrared limit. 

We proceed with the renormalization of the energy and the pressure associ- 
ated with the quantum field. 



Renormalization of the energy. 

krfn 



In the limit n 



El = 


/. 

27r2V 




o 

ry 




27r2A;' 






0. 




oo, we have 





+ (47) 

(48) 
(49) 



El = lim E^ = ^; (50) 
- (51) 



iJ^=lim_-^ = 0. (52) 
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Hence, the renormalized energy is: 



(53) 



Renormalization of the pressure: 



67r2 



- ei^i^.-^^; (54) 



o2 



(55) 



9n ^ Gir'^k' 

0. (56) 



9n 



In the Umit n — > oo, we have 



P0= limP," = g; (57) 



q2 pn 

P2 = lim = 0. (59) 

Hence, the renormahzed pressure is: 

P'r = |^,^[2k-l-,-2k + 1^^.0. (60) 

Again this is a result that we could expect from the begining since the de 
Sitter space-time has maximal symmetry. As consequence, there is a time-like 
Killing vector, and no particle can be created, and the energy and pressure 
associated to the quantum field must be zero. 



5 Inflation: the general case 

Let us consider now the general inflationary case. The scale factor reads, 

a = ao?7'/(i+3"^ (61) 
There is inflation if a < —1/3. The Klein- Gordon equation reads, 

0" + 2/3^+A:2</> = O, /3=:r^. (62) 
r] 1 + 3a 

The general solution is 

4> = ii-i^c^Hi^\kin)+C2Hl^\kii)^, ^=\-P- (63) 
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As before, we can choose ci in order to have the Bunch-Davies vacuum in the 
hmit k —^ oo, and C2 = 0. Hence, the final expression is 

cj.^T^'^cHi^Hkr,), (64) 

where ci does not depend on k. 

The corresponding expressions for the energy and pressure are 

p = A^'-^-'f' r k4Hl'\{kfj)H^^\{krj)+Hl^\krj)Hl^\krj)\dk, {^^^^ 



V-'^ r k^lHl'\{k^) Hi'\{k^) - iHl^Hkrj) ff(2)(fcr;)ldfc(66) 



Jo i -J 

where A = AttcI- 

The conservation law 

p' + 3^{p + p)=0 (67) 
V 

is verified. 

As in the preceding cases, there is no divergence in the infrared limit. Let 
us find first the divergences that may occur at the limit k ^ oo. Using the 
expansion (|17ll8p . we must go until order k^^ in the brackets of the integrands 
of (|65l66p : due to the term k'^ in the integral of the energy and pressure, all 
terms up to this order in the product of Hankel's functions are divergents. In 
particular, the term of order k~^ gives a logarithmic divergence, the term k~^ a 
quadratic divergence, and the term k^^ a quartic divergence, and all subsequent 
terms in the expansion are zero in the limit fc — > oo. Note that in both examples 
treated before, the logarithmic divergence is absent. 

Using (fTTlTSl) . we find 

TTX 

Q,{x)P,-i{x) - iQu{x)Qu-i{x)]] (68) 

h12^{x)hI^\x) - —[iP,-iix)P,{x) + P,.iix)Q,ix) ~ 
nx 

PAx)Q.^i{x) + iQ,(x)Q,_i(x)]; (69) 

H^^\{x)H(^\x) + hI^\x)hI^\{x) = ±[Q^(a:)P,_,(:,) - 

nx 

PAx)Q.-i{x)]. (70) 

Using the expansion for the functions P and Q, the expression for the energy 
and the pressure become: 

. s-fl 9 2 f„ (2iy - 1)2 
Ek = Arf''^^'^^^-x{ 2 ' ^ ' 



2 



y Ax 
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Pk = An 



,2i/-2/3-2 



2 (2 4jy2-12j/ + 5 



—x< - 

TT 1 3 



12x2 

(21/ - 3){2iy - l){2v - 9){2iy + 1) 
64x4 



+ • 



(72) 



where x = kr]. It is clear, at this approximation, that there are a quadratic, a 

quartic and a logarithmic divergences in the ultraviolet limit. Let us compute 
the different components to be subtracted from the above expressions, as we 
have done before. 

5.1 Computation of the energy 

The energy reads, 



where 7 = 2(1/ — /3). As consequence, we find 

£;« = r;f nfc2 1 hI^}-^ {nkrj) {nkrj) + Hl^\nkr]) 7J^2) (^^kr]) 



(73) 



(74) 

where X ^ X = nkrj. 

Consequently, the expression for the energy can be rewritten as: 

E]: = v'^-'kfix), fix) = xl^HiZ{x)Hi%ix) + if«(x)if(2)(^)|. (75) 
Hence, we find: 

El = \i^j^'^-^kx{H^^\ix)Hf\ix)+Hi^\x)H'^\x] 

Hi'Mx)Hl'},-HWix)H('\x) 



El = lim <! - '^%lrf-^x^ 

X->CX3 I 2k 



(76) 
(77) 



2v-\ 



Et = lim -—r^v-^-' xH i2v + l)/fWi(x)F(l\(x) + i2v - 3)hJ'\x)hJ'\x) 



x—>oc 4fc'^ 

- 2x fci) ix)Hl'}, (x) + Hi'}, (x)// (2) (^) 

Now we are read to compute E^, El, Ef and E).. Using first the relation 

Hi},ix)Hi'\{x) + Hi'\x)Hi'\x) = 

4i/4 - 4j/ + 1 481/"* - 96i/3 + 24i/2 + 24i/ - 9 



(78) 



-I 

TTX [ 



2 + 



+ 



4x2 64x4 

Hi'\x)Hi'},ix)+Hi},ix)Hi^\x) = 



(79) 
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2 r 2i/ - 1 384v^ - 576i/2 - 96^ + 144 
■ 4- 



TTX \ X 3(8a;)3 

(4i^2 - l)(4i/2 - + 3)(2t/ - l)(16z/'* - 32i/3 



1 



3(8x)5 
140z^^ + 456i^ + 1736) 
we obtain from (l76l). 



TT 



On the other hand, using 



Bi^ + 4 



3 (-64I/3 + 96i/2 + 16t/ - 24) ' 



2 64a;4 
and inserting in (j77[) we find 

In the same way it is possible to show that 



where 



/(z^) = ( 128j/^ - 576z^^ - 256t/^ + 5968i/^ + 2\21v'^ 



12288 

- 1020401/'* + 144321/^ + 554155;^^ - 1978j/ - 10794 



From (|17ll8p it can be written, 

4fc (2z/-l)2 48z/4 _ 96z/3 + 24J/2 + 24z/ - 9 ^ 
Ek = —rp + — TT-i + 57rT3 

Consequently 

^ ~ 2^'= " 349152^^ ' 

where 

h{v) = -256i/^° + 1280t/^ - 64;/* - 12192^^ - 3657i/^ 
- 225352;/^ - 571761/^^ - 1241334;/^ 
+ 5954975i/2 + 37044;/ - 14903. 
Finally, we find 

^ 491527r Jo k ' 

revealing that a logarithmic divergence remain. 
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5.2 Computation for the pressure 

The expression for the pressure is 

Pk = V'k' !^HlZikv)H!^Ukv) - lHi'\krj)Hi'\kv)^ . (90) 
From the definition (fT3|) it comes out, 

P[' - V^k'n S^Hi'\{x)Hi%{x) ~ \h'^J-\x)H'^:^\x)\^ , (91) 
where x — nkrj. Hence, 

= hm 7y^-ifcx(i?W (x)il(i\(x)-ii/«(a:)i?W(.T)l, (92) 

3Hi'-\{x)Hi'\ (x) + iJ (x)iJ (a;) 



Ax 



Hi'\x)Hl'\{x)+Hl'-\{x)Hi'Hx 

(12^2-3)i/«i(x)i?f_\(x) 



(93) 



+ x{-4v-10) 



Hl'\x)Hl'\ix) + Hl^\{x)Hi'Hx) 



8x' 



Hl'},{x)Hl'2,ix)-Hl'Hx)Hl^\x) 



(94) 



Using ini) and IHl) we find, 



2 r2 



i?i'_\(x)ilf_\(x) - \H('\x)Hi'Hx) 

3 TTX \ 3 

Av^ ~l2v + h IQv^ - 9Qv^ + 104J/2 + 2Av - 27 



12.t2 



leading to 



In a similar way, we find 
Pi 



64x4 



4!.2 - 12jy + 5 



Pi 



737287rfc3 



,7-5 



,7-5 



(95) 
(96) 

(97) 
(98) 



where 



g{i^) = 256i.^° - 512i/^ - 3456j.^ + 126721/^ + 2391681/^ - 642912;.^ 

- 2783281^4 - 67762i.3 - 246453i.2 - 150032i/ + 29838. (99) 
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From (pni we have, 

(100) 



Consequently, 



pren _ p _ pO _ pi _ }_ p2 _ K^) 7-5 /1 ni \ 

^ - 2^'= " 147456^fc3^ ' 



where 



Z(z^) = -256i/^° + 512i/^ + 3456j/^ - 12672;^^ - 239168i/^ + 642912:/^ 

- 204600;^'^ + 235256;^^ + 7256851/2 + 260624zy - 154254. (102) 

The renormaUzed pressure takes again the form, 
with also a logarithmic divergence. 



6 Renormalizing the general inflationary case 

The presence of the logarithmic divergence in the original expression for the 
energy and pressure seems to be a general feature of the problem of particle 
production in cosmology. See, for example, the reference [IJ. One important 
aspect concerning the logarithmic divergence is that the method of renormaliza- 
tion employed seems inefficient in the sense that not only it does not eliminate 
the logarithmic divergence in the ultraviolet limit, but also it adds a logarithmic 
divergence in the infrared limit. Hence, we propose an alternative scheme to 
cope with this indesirable feature. 

The proposed procedure is the following. First of all, we write the expressions 
for the energy and pressure as 

p - V"^"'''^'^"^'{i^^-\(a;)i^^i\(^)+ff^'Ha;)i^?^(a:)|rfa;, (104) 

p = V'^-2/3-5^°°a;4|i?W^(x)i/fJi(x)-ii/W(x)ij(2)(^)|d^^ 

where x = krj. This simple redefinition is well justified if 77 0, that is, the 
energy and pressure are computed out of the singularity. But the computation 
can be carried out as near of the singularity as we want. Hence, we have 

E, = A7j''^-'^-'^x'\Hl'\{x)Hl'\{x)+Hl^\x)Hi'\x)\, (106) 



P, ^ A^''^-'^-'x'lHi'\{x)Hl'\ix)~^H(J\x)Hi^ (107) 



14 



The divergencies can then be written as 



(2i.-3)(2i/-l)(2i.-9)(2i/ + l) ■ 



We subtract the (divergent) expressions for energy and pressure as 

Sr = Ek-E'^-El-Ei"^, (110) 

pren ^ _ pO _ ^2 _ plo9 ^ (^^l) 

with the definitions 

^ V^-^/3-5 6(2.-3)(2. 1)^(2.+ 1) 

TT 64 X'^ 

Pjo. ^ ^^2.-.,-5l (2^ - 3)(2:. -^1)^2. - 9)(2. + 1) _ ^^^3^ 

where an extra term has been added in order to assure that the divergence 
is eliminated in the ultraviolet limit without creating a new divergence in the 
infrared limit. The parameter a defines the efficiency of this mechanism, and 
the final result must be independent of it in order to guarantee the consistency 

of the procedure proposed. 

The energy and the pressure are now given by 



poo poo 
Jo Jl/a 

t-OC POO 

/ x\Pk-P^-Pl)dx- x^P]:Hx. (115) 

Jo J\/a 



^o^ 
ll/a 

Essentially we must compute an integral of the type 

-ax 



f ^—^ dx, (116) 

Jl/a X 



showing that the logarithmic divergence can be eliminate in the ultraviolet limit, 
when T 00, remaining a final result independent of a. In fact, this expression 
can be defined as 

1= / —dy, (117) 

Ji y 



where y = ax. The result is 

/ = Inr 



/ —dy. (118) 

Ji y 
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The first term eliminates the ultraviolet logarithmic divergence in the limit 
r ^ oo while the second one remains finite in the same limit. Moreover, the 
result is independent of cr. This procedure can be easily generalized if the cut-off 
is introduced as 

(l_e---")^, (119) 

where n is any positive number. 

As a result, we remain with the expressions for the energy and pressure 

p^^^ = Airf^-'^'^-^h = (120) 

p-" = A2?72--2/3-5j^ ^ ^^^^HOt^j^^ (^21) 

where Ai and A2 are constants and Ii and I2 are the (finite) integrals resulting 
from the renormalization procedure sketched before. 
Notice that, since the Ricci scalar can be written as 

R = (x rj-^^ , (122) 

the pressure and energy density are proportional to the square of the Ricci 
curvature. Hence, when the null energy condition is violated (a < —1), the 
energy of the created particles diverges as the singularity is approached, while 
in the case the strong energy condition is violated but the null energy condition 
is satisfied, the energy of the created particle is very high in the beginning 
decreasing as the the universe expands. 

Hence, using the renormalization procedure proposed here, the energy of 
created particle may be enough to avoid the big rip, in opposition to what has 
been stated in [3]. In fact the ratio of the renormalized energy of the created 
particle to the energy density of the dark energy fluid px is 

'- — cxr;-«— , (123) 

Px 

going to zero when —1/3 > a > —1 and to infinity when a < — 1. The case 
of the de Sitter space-time a = — 1, is just the limiting one, the energy of the 
created particle neither decreasing nor incresing, remaining always equal to zero. 



7 Conclusions 

In this work we have re-analyzed the problem of particle creation in a universe 
dominated by dark energy. Special attention has been given to the case where 
dark energy is represented by a phantom fluid. This problem has already been 
addressed in reference [3] . There the conclusion was that the energy associated 
to the created particles, as the universe approaches the big rip, goes to zero in 
spite of the fact that the number of created particles goes to infinity. 
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The conclusion presented in reference [3] has been obtained by introducing 
a cut-off in the energy integral at Planck's scale. This seems to be natural be- 
cause it can be expected that the usual Klein-Gordon equation must be modified 
beyond the Planck's scale, and some proposals in this sense suggest the intro- 
duction of an exponential decreasing term with the wavenumber fc, suppressing 
the contribution due to the transplanckian scales [IT)] . 

Here, instead of introducing an upper cut-off, we have chosen to renormalize 
the energy integral using the n wave method described in reference [S] . First we 
have shown that the employement of this method leads to the expected conclu- 
sions concerning the radiative and de Sitter universes: the energy expression is 
divergent, but after renormalization it becomes zero. 

Applying the method to the general dark energy model, a difficulty appears 
due to a logarithmic divergence in the ultraviolet limit. We have proposed a 
method to cope with this difficulty, introducing a new cut-off. The final result 
is independent of the cut-off. Moreover, the final expressions for energy and 
pressure are covariantly conserved. 

The result for the energy density differs from that found in reference [3], 
since now the energy density of the created particles diverges as the the big 
rip is approached, and becomes the dominant component. Hence, following 
this result, quantum effects can be effective to avoid the big rip. In fact, the 
renormalized energy comes out to be proportional to the square of the Ricci 
scalar. Another consequence of the computation made in the present work 
is that in the non phantom dark energy models the renormalized energy goes 
asymptotically to zero. The de Sitter case is the separatrix of these two different 
behaviour. 

An important point concerning this result is its unicity. A comparison must 
be made with other renormalization methods, for example those described in 
and references therein. We hope to address this problem in a future work. 
However, the absence of any final dependence on the cut-off employed to renor- 
malize the energy, together with the fact that the renormalized energy and 
pressure obeys a covariant conservation law, indicates that this result is quite 
consistent. 



Acknov^rledgement: We thank CNPq (Brazil), FAPES (Brazil) and the brazilian- 
frcnch scientific cooperation CAPES/COFECUB for partial financial support. 
E.G. A. and J.C.F. thank GreCO, lAP, France, for kind hospitality during part 
of elaboration of this work. We thank Olivier Piguet for the critical reading of 
the manuscript. 

References 

[1] N.D. Birrel and P.C.W. Davies, Quantum fields in curved space, Cam- 
bridge University Press, Cambridge (1982). 



17 



[2] A. A. Grib, S. G. Mamaycv, and V. M. Mostc- panenko, Vacuum Quan- 
tum Effects in Strong Fields, Friedmann Lab. Publ., St. Petersburg, 
(1994). 

[3] A.B. Batista, J.C. Fabris and S. Houndjo, Grav.&Cosmol.l4, 140 (2008). 
[4] D.N. Spcrgcl ct al, Astrophys. J. Suppl. 170, 377 (2007). 
[5] R. Caldwell, Phys. Lett. B545, 23 (2002). 

[6] J.D. Barrow, A.B. Batista, J.C. Fabris and S. Houndjo, Phys. Rev. D78, 
123508 (2008). 

[7] J.D. Barrow, Class. Quantum Gravity 21, L79 (2004). 

[8] Ya. B. Zeldovich and A. A. Starobinsky, Sov. Phys.JETP 34, 1159 (1972); 
M. Bordag, J. Lindig and V. M. Mostepanenko, Class. Quantum Grav. 15, 
581 (1998); Yu. V. Pavlov, Theor. Math. Phys. 138, 383 (2004). 

[9] L. Parker, Phys. Rev. Let. 21, 562 (1968). 

[10] M. Lemoine, J. Martin and J-Ph. Uzan, Phys. Rev. D67, 103520 (2003). 

[11] Yu. V. Pavlov, Grav.&Cosmol. 14, 314 (2008) 



18 



